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INTRODUCTION 

The linear programming model has been designed for the 
evaluation of technologies and policies and includes the full 
range of interfuel substitutability, including substitution 
between electric and non-electric energy forms. It encompasses 
the entire energy system including all resources and demand 
sectors. Since the range of interfuel substitutability that 
is feasible depends on the supply and utilization technologies 
that are available, the model is constructed around these 


technologies. The technology related parameters which appear 


explicitly in the model are the efficiencies of energy conversion, 


delivery, and utilization devices; the emissions or environ- 
mental effects produced by the devices; and their cost. 
The important characteristics of the linear programming 
model may be summarized as follows: 
1. The model ancompasses the entire energy system in- 
cluding all alternate resources and both electric 
and non-electric demande. 
The full feasible range of interfuel substitutability 
is reflected. 
Technical, economic, and environmental characteristics 
of energy conversion devices (both supply and utiliza- 


tion) are incorporated. 


The load-curation characteristics of electrical de- -> 


mands are included. (This is a very important element 
that is frequently ignored in energy systom modelling.) 
Supply, demand, and environmental constraints are 
specified exogenously. 

Optimal supply-demand configurations are indicated by 
the model. The optimization may be performed with 
respect to cost, resource consumption, or environmental 
effects in a given planning year, or with respect to 
some arbitrary combination of these factors. It should 
be noted that the optimization ia performed on the 
basis of annual cost in the planning year rather 

than a minimum present worth over some number of years 
or over a planning increment. Dynamic characteristics 


of the energy system may be introduced by applying 


the model sequentially to a series of development perioda, 


with the growth of demands and equipment turnover 
specified exogenously. 

The model may be applied to regional energy planning 
or at the national level using national average param- 
eters. The physical siting of power plants is not 
considered in the model, although population density 
in the vicinity of a class of power plants may be 


reflected by appropriate weighting of the environ- 


mental effects. 


The model is sufficiently general, that it may be applied 
to regional energy planning with imports and exports of fuels 
fixed as exogenous supply and demand categories, or to energy 
system planning at the national level. In some respects the 
analytical model can be quantified with more precise costs and 
environmental constraints at the regional level, although re- 
source constraints are more easily defined at the national level. 
The appropriate area to be addressed in a regional study depends 
on the planning objective. For air pollution studies the air 
quality region will be the area of interest while for other pur- 
poses a river basin, megalopolis, or city might be the preferable 


region. 


Both near- and long-term planning horizons may be addressed. 


In the near term, the model would be applied to indicate the 
optimal configuration of the new supply and utilization systems 
installed over, say, a 5-year period. The optimization would 
be performed with respect to either total cost or annual cost 
in any given year. The growth of demands and the turnover of 
existing equipment would be considered in developing the supply 
and demand constraints that would operate in the planning year. 
For long-range strategic planning, the model may be applied to 
analyze optimal configurations of the energy system in say, the 
year 2000, where it may be assumed that much of the existing 


System will be obsolete. 





For purposes of technology assessment it is useful to 
first establish & base case with the model in which the tech- 
noiogy is excluded. By admitting the technology in a sub- 
Sequent run the impact of that technology on resource consump- 
tion, environmental effects, and cost may be determined. Policy 
questions may be addressed in terms of their impact on avail- 
able technologies and price structure and on the supply, de- 
mand, and environmental constraints that operate on the energy 
system. 

DESCRIPTION OF THE ANALYTICAL MODEL 

The energy system, on a regional or national level, may 
be represented in & network formet as shown in Figure 1. The 
network in this case is quantified with the energy flows for 
the year 1969 from alternate resources through the various 
energy conversion and delivery activities to Specific end uses. 
Each link in the network represents a process or mix of processes 
used for a given activity such as the refining of crude oil. 
Costs and evironmental effects may be assigned to each link. 
Similar energy system diagrams have been developed for the pur- 
pose of technology assessment for selected future years, (1) 
Examination of the energy demand sectors at the right-hand side 
of the diagram indicates the degree of disaggregation that is 


required in the demand sectors. To evaluate the potential 
















for example, it is 





benefits of a hioh-efficiency air-conditioner, 








necessary to define the demand level for this specific end use 














and to estimate the degree of implementation of the improved sys- 






tem whether such implementation comes about purely from competi- 








tive market forces or from a policy or regulatory decision. It 






is also necessary, particularly in the example of the air- 










conditioner, to establish the electrical load-duration charac- 


















teristics of the demand sectors, as the type of generating 





equipment installed to serve that demand depends strongly on 





the load factor. 











It is possible to develop energy system diagrams by 
specifying the fuel mix and technology mix at some future ¡ int 


in time on the basis of available forecasts and trend analysis. 









In view of changing prices, demand levels, and fuei supplies 








it is desirable to have at an analytical model that will 







indicate optimal fue: mix configurations within the resource, 





environmental, and demand constraints that operate on the energy 








system. The linear programming model provides this function. 











Rather than using a network algorithm for further analysis and 





optimization, the energy system diagram may be transformed to 








a simpler structure. Figure l indicates that a given resource 









may be converted to electricity and one, or in some instances 









several, general-purpose fuels. Rather than to reflect these 





in a network structure it is convenient in the linear program 
to consider them as alternate supply categories subject, when 
appropriate, to a single resource constraint. 
The linear programming model is formulated about the 
classical transportation problem of determining the optimal 
routing of a product, in this case an intermediate energy form, 
from a set of n supply nodes to m demand modes where a cost and 
set of environmental impacts are identified for a unit of energy 
passing over each of then x m possible paths. The typical 
linear programming representation of the transportation problem 
is modified by the inclusion of efficiency coefficients in the 
Supply and demand constraints and is augmented by additional con- 
straint equations reflecting the environmental factors as well 

as certain technical features of the energy system. A graphical 
representation of the basic model and a definition of terms is 
given in Figure 2. 
The n supply nodes and m demand nodes can include all 
supply and demand categories. To provide a feasible path between 
a supply and demand catetory, both a supply and utilizing tech- 
nology must be identified. Fora given path, j, a resource Su 

is converted to intermediate energy form, X, at an efficiency, 
Bud: In turn, the intermediate energy form is used to satisfy 
demand D, at an efficiency Gu: A cost c, and set of environmental 
effects, included in f, are also defined per unit of inter- 


mediate energy form. 




























For a complete discussion of these techniques see Dantzig'^! 
(3) 







The mathematical formulation of the model is as follows: 








and Wagner. 
Minimize 














The 4 equations that /augment the supply and demand con- 
Straints include the environmental constraints and equationa 


that constrain certain energy flows to reflect technical fea- 











subject to 






tures of the energy system. The latter group of equations in- 









clude: 









l. Off-peak constraints that Specify the maximum amount 


of energy available from each central station electric 
















Source to serve off-peak electric or thermal demands; 


?. energy balance constraints that specify the amount of 


by-product energy available from Specific sources to 













serve other demands; 


3. storage balance eguations that ensure equality between 


the amount of energy supplied to storage and that de- 


x, 2 0. (1) 





















The summation in each constraint equation is over only 


those indices j that represent admissible intermediate energy 
livered from storage including losses; 


4. endogenous demand constraints by which portions of 


central station electric demands can be reassigned 


forms for the supply category, demand category, or other con- 
















straint. The full linear program array thus has k = n+m+ł con- 


















straints and these may be represented as: 






internally to categories with different load factors: 






and 


5. special constraints on individual intermediate energy 


forms. 


a x, sb i = lk, (2) 
















i T ve . int LO f la- 
where the aij include e. Be) and fu in the previous formu 


tion and the b. correspond to the S. D and Bas This problem 



















may be solved by the application of any of several algorithms . 
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illustrate an off-peak conatraint equation consider a 
tation supply category, such as a gas-turbine generator 
on fuel oil, which can furnísh electrical energy to 

1 categories, a peak electrical demand with a load factor 
d an off-peak electrical demand that can use energy at 
that it is available. If these demands were the only 
could be served by this supply cacegory, a peaking 

t equation would be required to ensure that the in- 

ower capacity to deliver an amount of energy Xi to the 

nd was not exceeded by the power involved in supplying 


€ energy x, to the off-peak demand. For every unit of 


2 
, up to that same amount can be supplied in the ideal 


2 and the off-peak constraint equation would be of the 


e > (3) 
x, + x, £0 


olution, x, were equai to X^ the installed power gener- 
acity would be base-loaded. 

general expression for an off-peak constraint equation 
ric supply category A where the off-peak demands have 
estriction and the supply plant can operate with a plant 
for the annual period is: 


x. $ O, (4) 


L 21) ) 


J 






a, LF for electric energy delivered to 
) j peak demands with a load factor 


T s PF, 











a,. = 1 for electric energy delivered to 
d off-peak demande, 














ay. = euj for thermal energy delivered to off- 
J peak demands. 






Additional correction factors may be applied to the coefficients 






to reflect other time restrictions on the use of off-peak power. 








If both seasonal and weekly peak and off-peak demands can be 






Satisfied by an electric supply category, then an off-peak con- 






straint equation is required for each period. 


An energy balance constraint is required when by-product 


heat is available from an electric generating device to serve 









other demands for thermal energy. If the conversion efficiency 








of such a device is e, the load factor of the thermal energy de- 






mand LF, and a fraction a of the waste power can be used ef- 






fectively, then the constraint is given by 








a (e-1) LF a. 4 









x. $0 (5) 









where Xi is the amount of electric energy produced and x, is 






the intermediate energy form, heat, supplied to another demand 







sector. 
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When provision is made for the storage of by-product or 
off-peak energy for use in other demand categories, the storage 
medium is defined as both a demand category that can accept that 
energy and a supply category that can, in turn, deliver the 
energy to other demand categories. A stcrage balance equation 
is required to ensure the proper baiance, allowing for any in- 
efficiencies that may be involved, between the energy supplied 
to storage and that supplied from it. 

The endogenous demand constrajints are developed in the 
following manner. Let demand category v, with a demand con- 
straint D, be served by two alternativo central station electric 
supply categories, which deliver intermediate energy forms 


x, and x and by a decentralized supply category delivering 


1 2' 
energy form x, Further assume that the primary load factor 
for D, is 0.5 but a fraction g of the intermediate energy form 
used at the primary load factor is used at a load factor of 
0.1. Also, let x4 and x, ropresent intermediate energy forms 
dzlivered at a load factor of 0.1 from two sources that could 


even be the same type of plant that delivers energy forms x, 


and vi Clearly, 


where E, is the amount of energy required by demand category v 


at the primary load factor of 0.5. The exogenous demand con- 


straint in this cas A = 
e, for d d 2" is 


d 
yy 0*9) x, + dv 0*9) x, + 443%, = D, 


anà the endogenous demand constraint is 


+ - — Lad 
gx, gx, X4 Xg O. (8) 


Any portion of 3 Let: i 
p à the demand H that is satisfied by the intermediate 


energy form x, does not involve any reassignment of energy to the 


endogenous demand category. This is because, unlike the case 


with central station electric Sources, no other equipment can 


be used to furnish the lew load factor portíon of the demand. 


The cost i 
st coefficients for variable X4 in this instance must be 


based on the composite load factor. 
Only one endogenous demand constraint is required for a 


given load factor category and portions of any other demands 
that have a component with that load factor may be assigned to 
this category in the same equation. 


The following example will illustrate the complete formula- 


tion of the model in a Simple case. 


Illustrative examples Consider three supply cate- 
gories (S,: ZŁ and $4). the first two of which are 


central station electric Sources, and two demand cate- 


gories (D, and D.). Let D, be a peak demand with a 


primary load factor of 0.5 but with a fraction g 
thet primary component required at a load factor 


0.1. Further, let D, be an off-peak demand with 


time restriction and assume that all ey; and dui are 


unity. Finally, let e be the cost and s, the air 


pollutant emission coefficients for variables s b 
Formulate the linear programming model for this prob- 
lem. 


Solution: The problem may be put in supply-demand 


matrix form as: 


Demands 


Di Endogenous 


(primary demand 


LF = 0,5) (LF = 0,1) 





Supply: central station,S, Xo 


central station,S, 


- 


decentralized,S, 


(The subscripts of variables in this example bear no 
relationship to those used in previous illustrations 


of constraint equations.) 








The solution to this problem is 


1 2 3 
x, + Se 
(1*g)x, + (1*g)x, 
*3 
“9X, + x, “GX, + Ks 
TX IX, + x, 
x, "BS, 


ew y | + Daf: + a,%, + a4x4 中 u 


21% + C7X, + €3*4 * "a". + C 5*6 


This example illustrates tl 
(eqs. 9 through 13), the endogenc 
the off-peak constraint (eqs. 15 
constraint (eq. 18), and the obje 
that variables x, and x, represer 
low load factor. In a larger prc 


could represent energy delivered 


Any electric gererating capacity 


side 


given by: 
£ 
$1 (9) 

* Ke * S, (10) 
x, + Se s S, (11) 

+ = 
x, D, (12) 
* Se + xa = D, (13) 
= 0 (14) 
< 0 (15) 
+ Se & 0 (16) 
(1499) / (1*9) x + x, 50 (17) 
+ Be Ze + a_x_ + Bee < B. (18) 
* Ge ie + CX, H c8xg = min z (19) 


ie supply and demand constraints 
us demand constraint (eq. 14), 
through 17), the environmental 
:ctive function (eq. 19). Note 
it energy delivered to Di at the 
tlem, each of these variables 
to à number of demznd sectors. 


installed to handle this load 


is also available for off-peak demands and this is reflected in 
the off-peak constraints for each central station supply cate- 
gory. 

If the value of B. the environmental constraint, is set 
high, then eq. 

The value may be set lower, however, to constrain the solu- 
tion within some total emission level. 

The coefficients in the objective function (eq. 19) would 
reflect the capital cost of the facilities used in the energy 
supply system as well as fuel and other operating costs. Social 
costs, or externalities, may also be included if these can be 
quantified. The capital and operating cost of utilizing de- 
vices such as home furnaces, air-conditioners, automotive engines, 
and blast furnaces may also be included. The typical cost co- 
efficient is thorefore made of the sum of at least two cost 
elements; a capital recovery cost, Ce: and an operating cost, 

ę* 

The capitai recovery cost, when the intermediate energy 
form is electricity, is calculated using the following relation- 
ship: 


_ X: (CRE) - 10° 
Ce ^ 3413 - 8760 - (PF) 


= capital recovery cost, $/10° Btu (electric), 
= capital cost of supply and/or utilizing facility, 
$/kW (electric), 

CRF = capital recovery factor, and 

PF = plant factor. 

The plant factor will generally be taken to be equal to 
the load factor so long as that provides for any plant shutdown 
period that may be required. However, if the load factor is 
1.0 or even near that value, the plant factor may be reduced be- 
low the load factor by whatever margin is desired. 

The operating cost includes the annual cost of fuels and 
resources consumed as well as the cost of other Operations such 
as fuel transport and distribution. In general, the fuel cost 
will include the capital recovery cost of any fuel extraction 
and treatment facilities that are required. 

It is important that the cost coefficients be defines in 
terms of the intermediate energy form represented by the variable 
that they operate on. When the intermediate energy form is 
electricity, for example, the capital recovery cost of central 
station electric plants wil! be defined properly using eq. (20), 
but the fuel costs must be multiplied by the inverse of the supply 
efficiency so that they refer to the electrical output rather 


than to the energy input to the plant. 
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The capital recovery cost for an electric supply category 
is à function of the plant factor which, in turn, depends on 
tho load factor of the demands that are served. When energy in 
delivered to a peaking or base-load demand category, the full 
capital cost of the generating capacity required to produce that 
energy is recovered and some additional capacity margin may even 
be included to account for plant shutdown periods that may be 
required. Thus, within the off-peak constraint that is specified 
for this generating capacity, off-peak energy may be delivered 
to other demand sectors that can utilize this energy form, with 
no capital recovery charge. The appropriate fuel cost is re- 
covered, however, in both peak and off-peak demands. The alloca- 


tion of costs in this manner does not imply any rate structure 


but merely ensures that the full costs of the energy System are re- 


covered. Thus, in the illustrative example, the cost coefficients 


€i: c: and c, would all include fuel costs. The coefficient c 


3 


would include no capital recovery cost element. Coefficients 
cy and c. would include the capital recovery cost, calculated 
for each case using the appropriate plant factor. 
In an optimal solution to the linear programming problem, 
each central station electric supply category is assigned to 
serve a specific demand or combination of demands. The electrical 


energy supplied from these sources, however, is a common product 





and the only significance to be attached to a specific assign- 
ment is the combined load factor of the demands that are served. 
Examination of the optimal solution will indicate which central 
station categories are base-loaded and which are employed for 
intermediate and peaking service. 

It is clear that there are no implicit assumptions built 
into the model regarding objectives or continued growth of energy 
demands. The analysis simply indicates the optimal configuration 
of the energy system within the constraints that are defined for 
a particular analysis. The annual cost of the overall energy 
system is minimized in each analysis. There will be cases where 
some individual demand sectors may bear a higher cost for ser- 
vice in an optimal solution than they would in some non-optimal 
assignment. The overall optimality condition, however, =nsures 
the existence of an internal rate structure that can redistri- 
bute costs such that some or all parties will be better off, and 
none worse off, than they would be with some non-optimal alloca- 
tion. 

The linear programming model has the capability of very 
high accuracy in determining an optimal solution; higher, perhaps, 
than is justified by the accuracy of the parameters that are 
defined. In particular, if the płanning horizon exceeds 15 years 


or so, the constraints and technological coefficients will be 


ole 





speculative. This indicates the need to recognize possible solu- 
tions that may be near optimal and to perform extensive sensi- 
tivity analysis in the planning process. In any event, the 
value of the particular analytical technique selected for this 
framework is not in its capability for high accuracy; but is in 
its ability to capture the essential structure of a complex sys- 
tem, such as that involved in the supply and utilization of 
energy, in a simple format where the assumptions and constraints 
are quite accessible. The admission of an objective function 
for optimization is another important factor in the selection 

of this approach. 

The detailed quantification and formulation of the model 
for a given planning task depends, of course, on the specific 
objectives of that task. It is important in any case to include 
all alternate energy resources and all demand sectors that 
compete for scarce resources at a reasonable level of detail; 
however, iräividual sectors, such as transportation or industry, 
can be included at a greater level of detail for studies directed 
specifically at these sectors. 

The level of disaggregation of supply and demand categories 
is optional and will, of course, determine the size of the linear 


programming problem. 


AN APPLICATION TO TECHNOLOGY ASSESSMENT 

In order to illustrate the general capability of the model, 
it has been quantified and applied to a long-range assessment of 
new energy technologies. A set of supply and demand constraints, 
technical coefficients, and costs were estimated for the year 
2000 and a series of optimization runs were made to evaluate the 
impact of a non-electric air-conditioning system and a battery- 
driven automobile by comparison with a base-case analysis in 
which these options were excluded. The data base and assumptions 
are described in full in the dissertation on which this paper is 


based (4) 


and will not be duplicated here. The results are pre- 
sented here primarily to illustrate the application technique and 
the output information that may be obtained. Additional runs 
using a range of technical parameters and cost estimates are 
required before any final conclusions may be drawn with respect 
to the technologies evaluated. 

The major air pollutants and high-level radioactive wastes 
were incorporated as environmental effects in this analysis. It 
was assumed that proposed EPA standards for central station power 
plants and automobiles would be achieved. 


The supply and demand categories identified for this study 


are shown in Table I, which also illustrates the optimal solu- 


tion for the base case in which the new technologies under consider- 


ation were excluded. Each feasible supply-demand combination in 
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Table I is represented by a variable in the linear programming 
model. The supply and demand constraints that applied througii- 
out the series of runs are given in units of WÉI Btu. Note 
that the runs and columns do not sum to the constraint vaiues 
since supply and utilization efficiency coefficients must be 
applied to the variables. 

The electrical system load-duration curve is divided into 
three blocks with load factors of 0.1, 0.5, and 1.0. The load- 
duration characteristics of each demand category, that apply 
in the event that they use central station electricity in an 
optimal solution are given in Table II. A primary load factor 
is defined for each alon; with the fraction required at the peak- 
ing load factor, 0.1. Data of this type is not readily avail- 
able and the paramcters indicated in Table II should be con- 
sidered as preliminary estimates. The demand constraint on the 
peak electric demand category is specified by an endogenous de- 
mand constraint that is a function of the particular demand 
categories served by electricity in an optimal solution. 

In the base case, the private ground transport category is 
served with gasoline and the air-conditioning demand is served 
by electricity. Some of the space heat demand is also served 
by electricity, thus filling in the winter valley. The water 
heating and desalination demand categories use off-peak energy 


but no other off-peak demands are available to be served. 
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Pumped storage is used to the limit of its constraint value for 


peaking service with gas turbines providing the balance of the 


peak demands. 


The pumped storage supply category has the highest marginal 


value in this base case. The marginal value represents the 


cost decrease that would be obtained for a unit increase in the 


constraint value for the supply category. The values are obtained 


from the solution to the dual problem. 

The assignmert of a central station Supply category to 
serve a given demand category is not of primary importance. The 
significant information is the combined load factor of all of the 
demands served by a specific electric supply category. In the 
base case both the Light Water Reactor (LWR) and the Liquid Metal 


Fast Breeder Reactor  (LMFBR) are base loaded as is all but a 


MY 
small portion of the coal-fired electric category. 


The asterisks i” Table I indicate variables that could 


enter the optimal basis at a cost reduction less than 10¢/10° Btu. 


Oil is in short supply in this run; however it has a relatively 
low marginal value indicating that alternatives are available at 
only a slightly higher cost. These alternatives are methane (ob- 
tained from coal) for air transport and coal for petrochemicals. 
Three additional runs were made in which the non-electric 


air-conditioner and the battery-driven automobile were admitted, 





as 





first indivic 
version effi: 
given in Tab; 
mitted with 
compared wit! 
non-electric 
current syst 
automobile w 
(including pi 
other componi: 
taken for a í 

When tl 
replaced the 
electric dem 
duction of a 
and radioact 

The in 
in a $43 bil 
system effic 
the emission 
carbon monox 
that are mor 


The addition 


iually and then together. The cost, overall con- 
"iency, emissions, and wastes for these cases are 
le III. The non-electric air-conditioner was ad- 

a coefficient of performance (c.o.p.) of 2.5 as 
h a value of 3.0 for the electric unit. The assumed 
c.o.p. represents a technological advance over 

ems which range between 1.0 and 1.8. The electric 
as assumed to have an overall conversion efficiency 
awer conditioning, batter charge and discharge, and 
snta) of 0.6 as compared with an efficiency of 0.2 
jasoline-fueled vehicle. 

he non-electric air-conditioner was introduced it 
electric device. The resultant decrease in peak 
ands produced a significant cost decrease and a re- 
11 emissions with the exception of carbon dioxide 
ive waste. 

troduction of the battery-driven automobile resulted 
lion decrease in annual costs. The overall energy 
iency increased by about 3.7 percentage points and 

S characteristic of the gasoline-fueled vehicle, 

ide and hydrocarbons, decreased while other emissions 
e characteristic of central station plants increasec. 


al electrical load in this case is handled by coal-fired 


atte 


4 * 


stations since nuclear Categories were already base loaded and 
employed at their constraint value in the base case. If thc 
additional loads were handled by nuclear capacity the only 
residuals increase would be of radioactive wastes. 

The cost reduction attributable to the electric automobile 
stems from the low fuei cost in central station power plants 
and the better utilization of capital facilities as compared 
with the base case. All of the central station plants are base 


loaded except for a small amount of pumped storage required for 


peaking service. 


In the final case, both technologies were admitted simul- 


taneously and the results are tabulated. The total resource 
consumption was lowest for this case, 205.7 x 101? Btu as com- 
pared with 220.7 x 1015 Btu in the base case. While the absolute 
values for energy consumption indicated in this series of runs 
may not represent an accurate forecast or projection, the re- 
lative changes attributable to technological shifts do provide 
valuable planning information. In general, it is felt that 
the assessment of a given techaology is not particularly sensi- 
tive to demand forecasting errors and thus accurate absolute 
forecasts are not essential. 

Throughout this series of analyses it was assumed that the 


new technologies could completely displace the systems appearing 


in the base case. It was further assumed that all battery 
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